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We investigate a new kind of dark solitons and vortices that can exist in the exciton-polariton 
condensates. These structures have discontinuity in the excitonic part of the polaritonic field and 
exist due to an intrinsic multiplicity of the solutions for the exciton density in the given optical field. 

Reported solutions are characterized by two very distinct localization scales, and hence are coined 
as dual-core dark solitons and vortices. 
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I. INTRODUCTION 

Photons in optical microcavities can strongly couple to 
excitons and form new quasi-particles: exciton-polaritons 
or simply polaritons These unique half-light half¬ 

matter states exhibit rich physical properties and attract 
great interest from condensed matter and optics com¬ 
munities. Inheriting a light effective mass from its pho¬ 
ton component, microcavity polaritons demonstrate high 
temperature Bose-Einstein condensation d, Q and su¬ 
perfluidity [1, At the same time, due to their exci¬ 
tonic component, microcavity polaritons interact much 
stronger than pure photons in typical photonic setups, 
and thus represent a competitive and promising plat¬ 
form for ultrafast signal processing applications. In¬ 
deed, a range of fundamental nonlinear effects has been 
investigated with microcavity polaritons, including low 
threshold bistabihty 0, i, polarization multi-stability 
and switc hing |^11| . parametric scattering and pattern 
formatio n [l3-ll^ . excitation of dark- |17l - l^ and brigth 
solitons [13117^ 

The existence of dark solitons is supported by the com¬ 
bination of the positive effective mass of polaritons with 
low momenta in the cavity plane, and the repulsive non¬ 
linear interaction of polaritons. The two factors jointly 
ensure stability of a high amplitude background, which 
forms the tails of a dark soliton [13]. In planar systems, 
such as polariton microcavities, these conditions also fa¬ 
vor formation of vortices, and the two ty pes of solutions- 
dark solitons and vortices, often coexist [^. In particu¬ 
lar, vortices can emerge as the result of a transverse in¬ 
stability development of a dark soliton [2l| . Spontaneous 
parametric scattering can also lead to formation of var¬ 
ious types of rotating and stationary vortex-antivortex 
lattice arrangements in microcavities [13, Isol - I^ . Exci¬ 
tation of vortices can be triggered by perturbations, such 
as an external pulse or polariton condensate scattering by 
an obstacle [^. The latter mechanism underlines sim¬ 
ilarities of polariton dynamics in microcavities with hy¬ 
drodynamics of super-fluids [s^j. The analogy has been 
further developed in the recent theoretical and experi¬ 
mental studies of hydrodynamical properties of vortices 


and dark solitons in polariton condensates De¬ 

spite the above analogies with hydrodynamics and op¬ 
tics, polariton solitons and vortices have an important 
distinctive feature: they are inherently two-component 
states, and the effective masses of the photonic and ex¬ 
citonic components differ by many orders of magnitude. 
The fact that there are two components is somewhat ob¬ 
scured, when polaritons are described by a single order 
parameter function, but in this work we demonstrate that 
there are easily accessible regimes and solutions, where 
this approach breaks down. 

In this paper, we demonstrate that the dual¬ 
component nature of micro-cavity polaritons leads to the 
existence of the novel type of dark solitons and vortices 
which have no analogues in hydrodynamics, nonlinear op¬ 
tics and dynamics of atomic Bose-Einstein condensates. 
Such solutions emerge as the result of the specific type of 
bistability of the system - so called intrinsic bistability. 
This phenomenon has been known for sometime, though 
very rarely investigated, and it refers to the situation 
when material polarisation is bistable in the presence of 
a constant applied field [doj] . In our system, the intrinsic 
bistability manifests in multiple solutions of the excitonic 
wave function for a given amplitude of the intracavity 
electric field. We show that the intrinsic bistability pro¬ 
vides the new and unique mechanism for localization of 
nonlinear excitations. We coin such solutions dual-core 
vortices and solitons. Recently it was noted that polari¬ 
ton vortices have two distinctive characteristic lengths of 
localization (healing lengths) in the photonic and exci¬ 
tonic components j41| . However, vortices and dark soli¬ 
tons reported in that work can be described within the 
framework of the polaritonic order parameter equation 
of the Gross-Pitaevskii type, and exist outside the pa¬ 
rameter range where the intrinsic bistability and our so¬ 
lutions are possible. We also note here the analogy with 
Abrikosov vortices in superconductors of the second type 
placed in magnetic field [13,113 • Abrikosov vortices have 
two very different scales: the so-called London length, 
which is the relatively large characteristic scale of the 
localization of magnetic field, and the coherence length, 
which is the small characteristic length of the variation 
of the superconducting electrons density. 
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II. MODEL EQUATIONS 

We adopt the well established dimensionless mean-field 
model describing dynamics of the exciton polariton con¬ 
densate in an optical cavity 0,0 in terms of coupled 
complex amplitudes E and "0 of the photon and exciton 
fields: 

dtE - iV\E + ( 7 i - i5)E = i'4, + Ej,- (1) 

— zcrV^')/'+ (72 + *1^ + (2) 

where time is measured in the units of the inverse Rabi 
frequency T = l/w^, scaling of spatial coordinates L is 
determined by the effective cavity photon mass rric'- L = 
7 i and 72 are the attenuation coefficients 
for the photon and exciton fields, respectively, 25 is the 
detuning between the exciton resonance and the cavity 
resonance, a is the coefficient describing the diffraction of 
the exciton field (relative to the diffraction of photons): 
a mclmex ^ 1 , is the exciton mass, Ep and qp 
are the amplitude and the detuning of the pump from the 
center of the gap between lower- and upper polaritonic 
branches. The coefficient of nonlinear interaction in the 
excitonic field is set to unity by the appropriate scaling 
of field amplitudes, see more details in Ref. (^ . 

For zero pump, Ep = 0, and neglecting dissipation and 
nonlinearity, the spectrum of low-amplitude zero momen¬ 
tum polaritons Ejip g-iujt jg gjygj^ | 3 y 

uj = ±y/lT¥. (3) 

This defines the gap between the lower and upper po¬ 
lariton branches, see Fi^T] Note, that Voronova et al 
assumed 5 = 0 in Ref. |4l| . which makes impossible ex¬ 
istence of the dual-core solitons discussed below. 


III. INTRINSIC BISTABILITY AND 
STATIONARY SOLUTIONS 

In this section we consider stationary solutions, E = 
ijj = in the conservative limit of negligible 

dissipation and zero pump: 71 = 72 = Flp = 0. Ampli¬ 
tudes A and dt solve the following set of equations: 

+ dy) A-\- {q + 5)A -|- 4 / = 0 , (4) 

(73^4'-b (g - 5 -b A = 0 . (5) 

Neglecting the diffraction term in the excitonic field, 
cr = 0, Eq. dS]) becomes the cubic algebraic equation: 


(g - (5)4'- 4'^-b A = 0 . ( 6 ) 

Generally, for a given amplitude A of the photonic com¬ 
ponent, it admits either one or three real solutions: 
4- = 4'™(A), TO = 1,2,3. 

We define coexistence of the three distinct states of 
exeitons for the same amplitude of the photon field as 


the intrinsie bistability of the system. We emphasize, 
that the intrinsic instability refers to the excitonic field 
only, and it is different from the conventional bistability 
of homogeneous solutions of the coupled system of Eqs. 
(HI) and (ID). 



FIG. 1: (Color online) Domains of existence of the non-trivial 
spatially homogeneous solution (shaded areas si, s 2 , /I, / 2 ), 
and intrinsic bistability domains (green areas si and s 2 ) in 
the plane of parameters (g, 5). 


The non-trivial spatially homogeneous solution of 
Eqs. (HIEI) with A{x,y) = Ah = const., 4'(a;,y) = 4'?i = 
const, is given by: 


4 - 


2 _ 

h 


q^ — 5^ — 1 


q + 5 


2 n 
" (g + (5)2- 


( 7 ) 


The domains of existence of this solution in the plane of 
parameters (g, (5) are indicated in Fig. [T] as four shaded 
areas sl,s 2 and /l,/ 2 . 

Analysis of Eq. shows that it admits three real so¬ 
lutions when 4'^ < (g — (5)/3. Combining this result with 
the solution for 4/^ in Eq. ©, we obtain the condition 
for the intrinsic bistability: 

2(g-<5)>^>0. (8) 

g -b 0 

In Eig. [T] the domains of the intrinsic bistability are 
marked as si and s 2 (green), while domains /I and /2 
(dark grey) correspond its absence. Homogeneous solu¬ 
tion Eq. 0 exists (4^^ > 0 ) across all these domains. 
Linear stability analysis shows that the homogeneous so¬ 
lution is dynamically unstable in /2 and s2 domains. We 
remark that it can be stabilized by including sufficiently 
high losses and pump, however such a strongly dissipa¬ 
tive case is out of the scope of this paper and will be 
discussed elsewhere. On the contrary, the homogeneous 
solution is linearly stable in si and /I domains. Hence 
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it can serve as the background for possible dark solitons 
and vortices. 

To describe localized solutions, we now consider spa¬ 
tially nonuniforin amplitudes A and 4'. For cr = 0, with 
the account of three possible solutions, d'm, of the al¬ 
gebraic equation Eqs. (HEl) can be written as the 
equation for the photonic component A: 

(d^ -l- dy^A -|- ((jf -|- (5)A -|- Tm(yl) = 0 , 771 = 1,2,3 . (9) 

Let us first discuss the ID case in which the fields are 
homogeneous along y-coordinate. Then, Eq. m is equiv¬ 
alent to the equation describing dynamics of a particle in 
a potential: 


dlA = - 


dFm 

dA 


( 10 ) 


where coordinate x plays the role of the effective ’’time”, 
and the potential is given by: 


FmiA) = 


3'StniA) {q-6)'i>l{A) , {q + 6)A^ 


■ ( 11 ) 


In the derivation of Eq. di) we used d/dA = d/d'Sm ■ 
{dA/d'i>rn)~^: where dA/d'^'m is obtained by direct dif¬ 
ferentiation of Eq. ®. 
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FIG. 2: (Color online) The effective potential Fm, Eq. (Hill , 
for q — 0, q = —1.8 and q = —2.2. The detuning is fixed to 
S = -2. 

In Fig. [2] the potential Fm{A) is plotted using different 
roots f['m(A) for a set of values of {q — 6). As follows 
from Eq. ®, the potential is multi-valued within the 
range of the photon field amplitude |A| <2[{q — (5)/3]^^^. 
As (g — S) tends to zero, the overlap area shrinks, and 
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FIG. 3: (Color online) The phase planes for g = 0, <5 = —2. 
Panels (a) and (c) illustrate the solution with two symmetric 
discontinuities at a; = ±1. A special case with the direct 
switch from the upper to the bottom phase plane is shown in 
panels (b) and (d). 


the potential becomes a single-valued function, see the 
bottom panel in Fig. 

The corresponding phase planes for g = 0 and S = —2 
[i.e. {q — S) = 2, see the top panel in Fig. [2] are plotted in 
Fig.lSKa). Each phase plane, corresponding to a separate 
root has one saddle fixed point. Crucially, for <7 = 0 
the equation for dt is algebraic, and hence the fk field can 
be discontinuous. At the same time, continuity of the A 
held and its derivative d^A is still preserved by virtue 
of the differential equation ([5]). In other words, within 
the range of amplitudes |A| < 2 [(g — (5)/3]^^^ where the 
three phase planes overlap, a phase trajectory can switch 
from one to another, provided that continuity of A 
and dxA is preserved. This allows us to construct var¬ 
ious localized solutions, for which the phase trajectory 
connects two saddle points from different phase planes. 

An example of one such phase trajectory is illustrated 
in Fig. Ha). At a; —oo the trajectory starts in the sad¬ 
dle point of the phase plane corresponding to the root 
'k 2 (A), see top panel in FigH^)- Within the allowed 
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range of A, the trajectory switches to the phase plane 
corresponding to the root 4'i(A) and then jumps to the 
4/3 one, where it terminates at the saddle point, see the 
bottom panel of Fig. [31Ja). The corresponding solution 
is shown in Fig. He) and it represents a dark soliton. 
The field A is continuous, but the field dt has discontinu¬ 
ities at two points, where the switches to different phase 
portraits occur. The phase trajectory can also switch 
directly between roots '^ 2 {A) and 4 ' 3 (A), as illustrated 
in Fig. [3])b). The corresponding soliton solution is plot¬ 
ted in Fig. [SDd). Comparing A and 4' field profiles in 
Fig. Eld), it is obvious that the soliton core in the A 
field is much larger than in the di* field, thereby forming 
a dual core dark soliton. We note, that the switching 
between the two roots occurs at a particular value of the 
field amplitude Aq, for which the continuity of d^A is 
preserved. For fixed values of parameters q and (5, this 
condition uniquely defines the structure of the soliton, 
in particular its localization lengths in each of the field 
components, and the value of 41 field at the discontinuity. 
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FIG. 4: (Color online) Transition from dual-core to conven¬ 
tional dark solitons. In panel (a) profiles of dark solitons are 
plotted for different values of {q — 5). Solid curves correspond 
to the case of fixed S = —2, dashed curves to the case of fixed 
q + S = —2. Panel (b) shows the dependencies of the inverse 
contrast of |4t| in a dark soliton on (q — S). The curve marked 
as I is calculated for fixed S = —2, the curves marked as 77 is 
calculated for q + 5 = —2, see the corresponding white arrows 
in Fig. [U 


It is instructive to consider how the above dual core 
dark solitons transform into conventional dark solitons, 
which exist for {q — 6) < 0, see domain /2 in Fig. [T] We 
note that 'I' 2 (^) = — 4 ' 3 (A), and therefore the profile of 
lib I for the dual core soliton in Fig. EJd) is continuous, 
see Fig.llja). In the soliton core, the field iTj reaches its 
minimum, however the minimum is always larger than 
zero and corresponds to | 4 > 2 (^o)| = |'I' 3 ( 7 lo)|- Decreas¬ 
ing (q — S) and approaching the boundary between the si 
and /I domains, the range of A, where the overlap be¬ 
tween the three phase planes occurs, gradually shrinks. 
As a result, the soliton background drops, as well as the 
minimum of iTj field in the soliton core. To characterize 
the discontinuity in the soliton core, in Fig. Hb) we plot 
the ratio between the minimal and the background values 
of 4> as functions of {q — S). Curves 7 and II correspond 
to two different paths in the parameter plane ( 9 ,( 5 ), as 


indicated with white arrows in Fig. [TJ Approaching the 
boundary between si and /I domains, this ratio gradu¬ 
ally reduces, but remains non-zero everywhere inside si 
domain, thus confirming that the soliton has a disconti¬ 
nuity of the 4' field in its core. 

Crossing the boundary between si and /I domains, 
the relationship between 4* and A becomes unique (no 
intrinsic bistability). Inside /I domain, the potential F 
and the corresponding phase plane now have three fixed 
points: the centre point at A = 0 and two saddles at 
A = ±|A^i|. In this domain, usual dark solitons cor¬ 
responding to heteroclinic orbits, exist. The profile of 41 
field becomes continuous and crosses zero in the soliton 
core. Therefore, the minimum value of |4I| becomes equal 
to zero for such solitons, and the ratio between the min¬ 
imum and background values is always zero, see Fig. 2) 



r r 


FIG. 5: (Golor online) Radial profiles of dual-core vortex soli¬ 
tons with m = 1 (solid curves) and m = 2 (dashed curves). 
In panel (a) vortex solitons with the discontinuity in 4* field 
in the centre (r = 0) are shown. The vortex solitons in panel 
(b) have discontinuity along the ring (r = 1). The parameters 
are cr = 0, g = 0, <5 = —2. 

The described above dual-core dark solitons with dis¬ 
continuities can be generalized to two dimensional vor¬ 
tices. Introducing polar coordinates (r, 0), and look¬ 
ing for radially symmetric solutions in the form A = 
/(r) exp(fm 0 ), d' = p{r) exp(im 0 ), where m is the topo¬ 
logical charge of the vortex, it is easy to write the ana¬ 
logue of Eq. ([9]) for the functions /(r) and p{r). Nu¬ 
merical solutions of this equation corresponding to the 
dual-core vortices (m = 1 ) localized at discontinuities in 
4' field are shown in Fig. El 

IV. STABILITY ANALYSIS AND DYNAMICAL 
EVOLUTION 

Structural stability of the dark soliton solutions with 
discontinuities on addition of the kinetic energy (disper¬ 
sion) term in the exciton equation (cr ^ 0 ) is an im¬ 
portant problem, which we address below. To investi¬ 
gate this issue we solved the ID version of Eqs. and 
([5|) numerically for small cr. The results are shown in 
Fig. El One can see that, as expected, the finite exciton 
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FIG. 6: (Color online) Dual-core dark solitons with non-zero 
diffraction in the excitonic component, a > 0. Panels (a), 
(b) and (c) show the profiles of A (dashed curves) and ^ 
(solid curves) fields for a — 0.25, a — 0.025 and a = 0.0025, 
respectively. Other parameters are: S — 1.55, q = —0.45. In 
panel (d) the instability growth rates are plotted as functions 
of q for the soliton with one discontinuity (solid curves) and 
a = 0.25, 0.025, 0.0025, and a soliton with two discontinuities 
separated by a distance Aa; = 1, cr = 0 (dashed curve). 


dispersion removes the discontinuity in the dl field, while 
the overall structure of solution is retained. Thus dark 
solitons with discontinuities are structurally stable and 
transform into dual-core solitons for finite a. To observe 
dark solitons and vortices in real physical experiments, 
the solutions must also be dynamically stable. We have 
studied the stability of the solitons by solving the corre¬ 
sponding spectral problem and by direct numerical simu¬ 
lation of Eqs. 0-0 with the noise added to the soliton 
profiles. The dependencies of the instability growth rates 
on the soliton parameter q are shown in panel (d) of Fig. [6] 
for different values of a. 

The linear stability analysis of conventional dark soli¬ 
tons reveals that they become unstable with respect to 
drift instability when q is increasing and the same insta¬ 
bility is inherited by the dark solitons with discontinu¬ 
ities. Increasing dispersion in the exciton model tends to 
suppress this instability, see Fig. [6Kd). 

Despite the instability of the dark soliton and vortices 
with very small diffraction in the exciton field, we demon¬ 
strate that similar structures can be experimentally ob¬ 
served in the decaying condensate. In the experiments 
polariton systems are dissipative having life time of or¬ 
der of 10 ps. We performed numerical simulations of the 
2D polariton system excited by a vortex laser beam and 
observed dual-core vortices in the decaying condensate. 
In this case the instability is suppressed by the dissipa¬ 
tion and the observation time is limited by the lifetime 
of polaritons. 

To study formation of vortices in the decaying conden¬ 
sate we performed a numerical experiment exciting the 
system by an optical vortex beam with the topological 


(a) (b) 



FIG. 7: (Color online) The variation of the pump amplitude 
Epo{t) in the numerical experiment is shown in panel (a). The 
parameters are qp = 0, 5 = —2. The inset shows the spatial 
distribution of the intensity of the pump. The bifurcation 
diagram for the spatially uniform backgrounds is shown in 
panel (b). 


charge one: 

Ep = Ep^it^ X 

X + iy 
2r ’ 

where r = + for our numerical simulations we 

chose ro = 1, Wp = 0.25. The corresponding spatial 
profile of the pump intensity is shown in the inset of 
Fig.m 


Wr, 


,tanh(!:^'l-tanh^" + "° 


Wr, 



FIG. 8: (Color online) Fields evolution in the numerical ex¬ 
periment: actual fields E and 4^^ panels (a) and (b), and re¬ 
normalized fields e = |i5|/max(|i?| and 4> = |^'|/max(|4'|), 
panels (c) and (d). Parameter values are: cr = 0, = 0, 
5 = —2. The inset in panel (b) shows profiles of 14^1 field at 
t — to = 6 and t — to = 15, the curves are marked correspond¬ 
ingly. 

The amplitude Epo of the pump was adiabatically var¬ 
ied in time, as shown in panel (a) of Fig. [T] We switched 
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the pump on at t = 0 and choose the initial ampli¬ 
tude large enough to bring the vortex background to the 
high amplitude state, see panel (b) of Fig. [7] showing the 
bifurcation diagram for spatially uniform states. Then 
we waited to obtain a stationary state, and at t = 100 
we started to decrease adiabatically the pump until at 
to = 400 it reached the value Epo = 0.43 close to the 
lower folding point of bifurcation diagram for the spa¬ 
tially uniform solutions. 

Then we switched the pump off, which has led to the 
excitation of waves that affected mainly the E filed. Due 
to dissipation the fields started to decay, see panels (a) 
and (b) of Fig. [S] In panels (c) and (d) of Fig.[8]the corre¬ 
sponding re-normalized fields are plotted. The structure 
of the vortex is clearly preserved as fields decay. 



FIG. 9: (Color online) Field profiles in the numerical experi¬ 
ment for the times t = 6 (a) and f = 15 (b). Dashed curves 
show the stationary vortex solution for the conservative prob¬ 
lem. 

One can notice that initially the core in 'k field is much 
narrower compared to the one in E field, as it should be 
in a dual-core vortex. With the decay of the background 
the both cores of the vortex broaden, but the 4'-core 
remains much narrower, see the normalized field distri¬ 
butions in panels (c) and (d) of Fig. [S] The onset of a 
drastic broadening of the 4'-core happens at the moment 
when adiabatic decrease of the vortex g-parameter trans¬ 


forms the dual core vortex into a usual one. In Fig.[9]the 
structure of a vortex in the decaying condensate is com¬ 
pared against the numerically exact vortex solution, both 
are clearly very similar. 


V. CONCLUSION 

We have demonstrated that intrinsic bistability of ex- 
citons leads to the existence of a new type of dark soli- 
tons and vortices. Without dispersion of the exciton field 
these structures have discontinuities in the exciton field. 
Arbitrary dispersion in the excitonic field removes these 
discontinuities and results in dark solitons and vortices, 
which have very different core sizes in the optical and 
excitonic fields. These double core vortices and solitons 
exist in the different parameter range relative to the pre¬ 
viously studied polariton vortices, which can be well de¬ 
scribed by the single polaritonic amplitude (an order pa¬ 
rameter) and have comparable core sizes in the excitonic 
and optical components [4lj . Formation of the dual core 
vortices in the decaying condensate under the realistic 
excitation conditions takes place despite possible pres¬ 
ence of the dynamical instabilities, since the characteris¬ 
tic time of the instability development is larger than the 
life time of the polariton condensate. 
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